Double-Exchange Model: Phase Diagram at Zero-Temperature 



David Pekker, 1 Swagatam Mukhopadhyay, 1 Nandini Trivedi, 2 and Paul M. Goldbart 1 

1 Department of Physics, University of Illinois, 
1110 West Green Street, Urbana, Illinois 61801-3080, U.S.A. 

2 Department of Physics, The Ohio State University, 
191 West Woodruff Avenue, Columbus, Ohio 43210, U.S.A. 

The analytical zero-temperature phase diagram of the double exchange model for classical back- 
ground spins as a function of the carrier density and Hund's coupling in the entire range of these 
parameters is presented. By constructing a continuum field theory we explore the possibility of a 
continuous phase transition from ferromagnetic state to a gently varying textured state. We find 
such a transition in and below two dimensions and show that the emerging stable state is a spin- 
spiral which survives the tendency towards phase separation into commonly considered phases, and 
is also energetically favored to the canted state, for low carrier density. 
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I. INTRODUCTION 

The Double- Exchange Model (DEM) describes the mo- 
tion of noninteracting itinerant electrons through a lat- 
tice of classical spins to which the electron spins are cou- 
pled. This model is relevant, e.g., in the area of dilute 
magnetic semiconductors, such a GaAs doped with Mn, 
which are important for spintronics applications pj. In 
these materials, the 5 = 5/2 local moments on the Mn 
sites are exchange-coupled to the hole carriers generated 
by substituting Mn for Ga. While a more realistic model 
would include the effects of disorder 0, arising from 
the random substitution of the Mn atoms, as well as spin- 
orbit coupling effects Q , a careful study of the DEM is a 
necessary first step. Another class of materials for which 
the DEM is relevant are the manganites, which show 
colossal magnetoresistance. Although electron-phonon 
interactions play an important role in the manganites 
through a Jahn- Teller effect j5j , some of the key features 
of the strongly coupled spin, charge and lattice degrees of 
freedom are captured by a DEM with lattice-distortion 
effects. 

Thus, the DEM is a paradigm for a wide class of ma- 
terials that show a strong coupling between the charge 
and spin degrees of freedom. Pioneering work on this 
model focused on the ferromagnetic and canted antifer- 
romagnetic phases at nonzero temperature 6] . Over the 
past few years, the DEM has received renewed attention, 
stimulated in part by the numerical studies reported in 
Refs. [tIIUj where the zero-temperature phase diagram of 
the simple DEM has been explored. More elaborate ex- 
tensions of the DEM — intended to move it closer to real 
systems by augmenting it with physical processes such 
as the super-exchange and Coulomb interactions, disor- 
der, and Jahn- Teller distortions, etc. — have also been in- 
vestigated over the past few years. For example, the 
stability of a spin spiral state in one extension of the 
model was examined in Refs. |9| and was contrasted with 
the stability of the canted state. The ground state of 
a model augmented with Coulomb interactions and large 



Hund coupling was addressed in Ref. . The instability 
of the homogeneous canted state with respect to phase 
separation, for large Hund coupling in a model that in- 
cludes tunable super-exchange interaction, was studied 
in Ref. The phase diagram of a three-dimensional 

model in the infinitely large Hund-coupling limit was 
studied in Ref. [12|, and near the Curie temperature an 
instability with respect to spontaneous translational sym- 
metry breaking was proposed, as was the possibility of 
phase separation. Recently, Ref. |13| addressed the sta- 
bility of the spin spiral state as the ground state of the 
DEM in the large- S limit. 

In spite of the considerable amount of theoretical and 
numerical work on the DEM, the continuous phase tran- 
sition from the ferromagnetic phase to a spin-textured 
phase remains incompletely understood. Considerable 
attention has been paid to the model in the infinite Hund- 
coupling regime, motivated by the fact that in the man- 
ganites, the Hund coupling is large, compared with the 
electronic bandwidth. The zero-temperature phase di- 
agram has been recently studied in Ref. [l4| within dy- 
namical mean field theory (DMFT) for the entire range of 
carrier concentrations and Hund coupling. It was shown 
that the stable ground state in different regions of param- 
eter space is either a ferromagnet, or a commensurate an- 
tiferromagnet, or some incommensurate phase with an in- 
termediate wavevector. Moreover, a second-order phase 
transition (from the ferromagnetic to the incommensu- 
rate phase) and a first-order transition (from the antifer- 
romagnetic phase to a region of phase separation) were 
identified. 

In this Paper, we address the magnetic ordering of the 
single-band DEM with classical background spins at zero 
temperature (ignoring orbital or charge ordering). In 
contrast to the work of Chattopadhyay et al. [l4[, we 
use a continuum field theory and a gradient expansion 
to determine the critical line of continuous phase transi- 
tions separating the ferromagnetic and textured phases 
in the parameter space of electron density and Hund 
coupling. This approach allows us to consider all pos- 
sible long wavelength textures that can emerge from the 



2 



DEM. We argue that the spin spiral is the energetically 
favored one from amongst this class of textures. As is 
well known, the ground state of the DEM is phase sep- 
arated at large values of the Hund coupling. However, 
by explicitly comparing the energetics of phase separa- 
tion for commonly considered textures, we argue that 
there is a region of phase-space in which the continuous 
transition survives the tendency towards phase separa- 
tion. We also show that the spin spiral state is favored, 
energetically, over the canted state, thereby making pre- 
cise the nature of the emergent 'incommensurate state' 
found in phase diagrams presented in earlier work 0, Il4j . 
We present the phase diagram of the model, after in- 
cluding the commonly considered candidate phases for a 
pre-emptive phase separation at high Hund coupling. In 
systems for which the Hund coupling is small (e.g., the 
cobaltates, diluted magnetic semiconductors, etc.) and 
the disorder is low, the approach used here may be of 
relevance. 

The present paper is organized as follows. In Sec. [H] 
a continuum version of the DEM is derived. The sym- 
metries of this continuum version are considered and the 
electronic degrees of freedom are integrated out, yield- 
ing an effective Hamiltonian for the background spins 
within a gradient expansion. A line of continuous phase 
transitions is determined, and it is shown that the spin- 
spiral state is the emergent stable state. In Sec. IIHI the 
issue of whether the phase transition from the ferromag- 
netic states to the spin spiral state is pre-empted by 
another transition (to either an antiferro/ferrmagnetic 
microphase-separated state or to a canted ferromagnetic 
state) is considered. It is shown that there is a region 
of the phase diagram in which the spin spiral state is a 
stable state. 



where the first, second and third terms respectively de- 
scribe electronic hopping, double-exchange and super- 
exchange couplings. Here, S» is the background spin at 
lattice site i, which we approximate as a classical vector, 
Jh/S* is the strength of the ferromagnetic Hund coupling, 
and Jaf / S 2 is the strength of the antiferromagnetic cou- 
pling between the localized spins. We shall explore the 
Jh vs. electron density phase-diagram, considering Jaf 
to be vanishingly small. However, for technical reasons to 
be explained later, we shall need Jaf to be nonzero (see 
Appendix . Earlier work has mostly focused on the 
Jaf vs. electron density phase-diagram of the model in 
the regime in which Jh is very large, so that the electrons 
are aligned with the background spins. 

It is convenient to transform to a new spin basis at 
each lattice site, so that the local part of the Hamilto- 
nian (i.e. the part that dominates at large Jh) is diagonal. 
To this end, we rewrite Eq. (|2.1|) in terms of the new op- 
erators {dj,^} and {ej^} which are, respectively, the 
creation and annihilation operators associated with the 
spin basis aligned and anti-aligned with the background 
spin direction at site j: 
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where the spinors fj and are defined in terms of the 
the polar and azimuthal angles of the background spin 
direction at site j, i.e., 8t and (pi, and are given by 



II. DOUBLE-EXCHANGE MODEL; 
ANALYTICAL STRATEGY 

The Double-Exchange Model describes noninteracting 
itinerant electrons moving on a lattice of static "back- 
ground" spins whose moments are typically large com- 
pared to that of the electron spins, and hence may be 
treated classically. The Hamiltonian is given by 



H = -t 



X] {°ia C ja 



h.C. 



~V S J" °ja a ^C ]h + -|f Sl ' S J ' t 2 - 1 ) 



j ab 



•:.i> 



7, 



7! = e-** 



_ e -i4>j/2 gin ^ 
e +i^/2 cos ^ 



(2.3a) 
(2.3b) 



7 and 7^ form an orthonormal local basis corresponding 
to the aligned and anti-aligned spin states. This map- 
ping is defined up to a phase factor \j , which is a gauge 
freedom In the continuum limit, the Hamiltonian is 
(up to the corresponding gauge transformation) given by 



H = -Jdx ip\x) i (d a - iA a f + i (dph ■ dpii) - J H | ^(x) - <^ f (x) j(<9 Q + iA a ) 2 <p(x) + ~ (dpn ■ dpii) + J H | y(x) 



+{V f (x) A Q (x) a Q ^(x) + ^ f (x) (9 a A a (x)) p(x) + h.c.} - J AF (d p h ■ fyn) 



(2.4) 



3 



where ^>(x) and <p(x) are, respectively, the field (anni- 
hilation) operators describing locally aligned and anti- 
aligned electrons and x is the position vector, n(x) is 
the unit vector along the background spin direction, and 
A and A are vector potentials originating from a Berry 
phase and are defined by 

A a (x) = icos(9(x) <9 Q 0(x) (2.5a) 

A a (x) = -d Q 0(x) + isin(9(x) d a <j>{x). (2.5b) 

Note that A* A a — dpii ■ dpn. Also note that in pass- 
ing to the continuum limit we have retained only the 
lowest-order terms in the gradient expansion of the field 
operators. As a result, the electron kinetic energy, and 
hence the band structure that corresponds to it, is of a 
simple parabolic form. 

A. Symmetries of the continuum Hamiltonian 

1. Local gauge invariance 

The mapping from a vector (representing a classical 
spin) in three-dimensional space to a spinor in SU (2) 
is defined up to an angle y(x), an overall phase factor, 
which is a U(l) gauge freedom i.e., under the simul- 
taneous transformations 

A(x) -> A(x) + Vv(x) A(x) -► e 2i ^ x ~> A(x) 

(2.6) 

t/»(x) -> e lx ( x ^(x) 0(x) -> e"^(x)0( x ) 

the Hamiltonian is invariant. 

2. Global spin rotation invariance 

The Hamiltonian is also invariant under rotation of all 
the background spins by a global angle. One can show 
that the effect of a global rotation is identical to a gauge 
transformation. For example, let us consider rotations 
about a;-axis by a small angle u>. The unit vector n(x) 
transforms to rt'(x) where 9'(x) ss 0(x) — cjsin<^>(x) and 
</>'(x) « </>(x) — uj cot 9(x) cos </>(x). Notice that the above 
transformation leaves dpfi-d^fi invariant, but the vectors 
A(x) and A(x) transform, up to first order in u, in the 
following manner: 

A a - Av-ZeJ^), (2.7a) 
2 \ sm J 

A a -> e - iucos ^ sine A a . (2.7b) 

On identifying the factor — u> cos 4>/2 sin 8 with the gauge 
parameter x [see Eq. (|2.6|l ]. we see that Eq. (|2.7a|l de- 
scribes a gauge transformation corresponding to rota- 
tions. A similar analysis can be carried out for ro- 
tations about the y-axis; for rotations about the z- 
axis the transformation is trivial. Therefore, global 
rotational-invariance corresponds to the gauge freedom 
in the model. 



B. Effective Hamiltonian 

In this section we derive an effective Hamiltonian H e g 
governing the spatially-dependent background spin ori- 
entation (texture), in the limit that the texture varies 
on length-scales much bigger than the inverse Fermi 
wavevector. To do this, we integrate out the electronic 
degrees of freedom, assuming that the zeroth-level de- 
scription corresponds to a system in the presence of a 
spatially uniform texture. The contribution due to any 
inhomogeneity of the texture is then treated as a pertur- 
bation via a gradient expansion. The demerit of this con- 
tinuum approach is that background spin configurations 
that change abruptly from one site to another (e.g. canted 
states or antifcrromagnetic state) are excluded from con- 
sideration. 

The effective Hamiltonian that results from this ap- 
proach is a functional of A(x), A(x) and n(x) and their 
derivatives. Working at fixed chemical potential /i, the 
effective Hamiltonian is defined as 

e -/3fleff[A,A,n, M ] =Jx>ipVip e -/3(#"[V',<P,A,A,n]- A1 A0 

= / VipVip e ff o[^,rf+i?d'/w,AA,fi] ) (2.8) 

where N = J dx (tin (x)^(x) + <p^ (x)y>(x)) , Hq is the free 
Fermi gas Hamiltonian, and the perturbation Hi is a 
functional of A, A and <9n, each of which has one spa- 
tial derivative [see Eq. H2.5al2.5bl) ]. and therefore is small 
in the sense of our approximation scheme. The form of 
-ffeff is constrained by gauge invariance. Thus, keeping 
allowed terms to quartic order in gradients, but for now 
setting Jaf to zero, we find the following form, arranged 
in increasing order: 

H cS [A(x),h(x),A(x)} 

— J G?xja(/i) d a h ■ d a h + b(fi) (<9 Q n • <9 Q n) 2 

+ c(p)F a/3 F af3 + d{p)\D a A a \ 2 + e((i)\D a Ap\ 2 } 
+ •••, (2.9) 

where the coefficients a(/x), b(p), c(/i), d(p) and e(p) 
are evaluated by computing the corresponding Feyn- 
man diagrams. Not surprisingly, the results are com- 
pactly expressed in terms of the following quantities: 
F a p = d a Ap(x) — dpA a (x) and D a = d a — 2iA a ; note 
that the combination D a A a is gauge invariant. 

The Feynman diagrams contributing to the lowest- 
order terms in the gradient expansion are shown in Fig.^ 
The amplitudes of these diagrams are both proportional 
to (d a h) 2 . Therefore, their coefficients, p\ and p2, add to 
give the stiffness of the ferromagnetic state p; when p goes 
negative, the ferromagnetic state becomes linearly unsta- 
ble. The dependence of the corresponding limit of stabil- 
ity on d, p and Jh is among the central results of this pa- 
per. The contribution from these two diagrams compete 
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FIG. 1: Lowest-order Feynman diagrams (1 and 2): solid 
lines represent aligned electron propagators, dashed lines rep- 
resent anti-aligned electron propagators, curly lines represent 
A a , and dotted lines represent d a ii ■ d a h. 
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FIG. 2: Ferromagnetic stiffness in dimension d for d = 1.5 
(dotted), d = 2 (solid), and d — 3 (dashed). Here, Jh = 1 
and Jaf = 0. Note that for d = 2, the stiffness becomes zero 
at \i = Jh but does not become negative. 



with one another. The first, gives a positive contribution 
to the energy, because a spatial variation of the back- 
ground spins decreases the hopping amplitude, via the 
Anderson-Hasegawa mechanism: t — ► tcos(©y-/2) where 
Qij is the angle between the nearest-neighbor spins i and 
j in the discrete version of this model. The second di- 
agram gives a negative contribution; spatial variations 
in the background spin orientation allow for mixing of 
aligned and anti-aligned bands, thereby lowering the en- 
ergy. The contributions to the stiffness [i.e. the coefficient 



of (d a n) ] are: 



Pi 



Pi 



2 1 + d TT d / 2 dT[d/2} 

+ Jn) 1+d/2 - (p - J H ) 1+d/2 e\fL - J H ] 
2 d + 1 n d / 2 (2 + d)J H r[l + d/2] 



(2.10a) 



, (2.10b) 



where T[-) is the Gamma function, ©[■] is the Heaviside 
step function, and d is the dimension of space. 

By examining the stiffness as a function of /x and d (see 
Fig. |2J) we observe that d = 2 is a threshold dimension, in 
the sense that the instability occurs for dimensions less 
than two (but not for dimensions greater than two). We 
emphasize that the precise location of the transition, as 
well as the threshold dimension, depends on the form of 
the bare electronic band structure, which we have taken 
to be parabolic. Corrections to the parabolic electronic 
dispersion relation would alter both the location of the 
transition and the threshold dimension. 

In particular, in two dimensions the contributions com- 



T 2 



32ttJ h 



(2.11) 



This shows that in two dimensions and at zero tempera- 
ture there is a critical chemical potential /i c = Jh above 
which the ferromagnetic phase loses stability and, as we 
shall see, undergoes a transition to a 'textured' phase. 
This critical chemical potential coincides with the bot- 
tom of the anti-aligned electron band. To investigate 
how the instability of the ferromagnetic state is resolved, 
it would therefore be necessary to raise the chemical po- 
tential above its critical value, which would begin popu- 
lating this band. However, as the gradient expansion is 
an expansion powers of qtexture/fcp Ppcr , it would not con- 
verge, as fcp Ppcr would be very small. We approach this 
dilemma by noting that these results were obtained in 
the absence of an antifcrromagnetic term in the original 
Hamiltonian 12.11 i.e., for Jaf = 0. The precise location 
t^c(J}i,d) of the instability of the ferromagnetic state is 
perturbed, and in general shifted to lower value, in the 
presence of a positive Jaf, thus creating a region in the 
(/i, Jaf) plane in which the ferromagnetic state has be- 
come unstable and yet \x is still smaller than Jaf, sothat 
the upper band remains unoccupied. This scheme opens 
up a region of the phase diagram in which our gradient 
expansion remains valid and, at the same time, a tex- 
tured state is prefered. 

In order to investigate the form of the (stable) textured 
state that replaces the (unstable) ferromagnetic state at 
chemical potentials immediately greater than the criti- 
cal one (given in two dimensions by [i = Jh), we ex- 
pand the effective Hamiltonian density to fourth order 
in gradients of n, and minimize it with respect to all 
textures that vary only on lengthscales longer than the 
Fermi wavelength. Via the extension to quartic order 
of the diagrammatic expansion described in the present 
section, and for the case of two dimensions, we find the 
effective Hamiltonian to be 



H c ft — 
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J H 
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32ttJ h 
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48ttJ3 

(J h + a*) 2 (2Jh - m) 
192ttJ3 



M-/ h + m 2 ) 

256ttJ£ 



(d a h) 



-F af3 F, 



a 6 



D a A \ 2 



J A-^\D a /x a \ 2 - J AF (d a ii) 2 , (2.12) 



96ttJ3 

where the terms associated with Jaf arise from the an- 
tiferromagnetic term in Eq. (|2.1(l . 

The details of making this extension to quartic or- 
der are given in Appendix [B] For fi larger than its 
critical value, the coefficient of the first term in H e ft 
(e.g. the ferromagnetic stifness) is negative, and there- 
fore it is favorable for the ground state to have a non- 
uniform texture. We now show that the remaining terms, 



which are fourth order in gradients and serve to restabi- 
lize the textured state, are all positive-definite whenever 
< fi < Jh- The coefficients of (<9 Q n) 4 and F a pF a p are 
positive for < fi < Jh, ensuring that these terms are 
indeed positive-definite. If we neglect the surface terms, 
the forth and fifth terms can be recast in the following 
form: 



fj, (Jg - fi 



64tt 



48ttJ3 



(2.13) 

In this form, the coefficients of each of these terms is pos- 
itive for < /i < Jh, and therefore all the fourth-order 
terms are indeed positive-definite. Keeping only the first 
two terms, it is easy to check from the differential equa- 
tion for the ground state, which follows from varying H e s 
with respect to the fields 9 and <fr, that a spin spiral state 
[e.g. 9 = ir/2 and <f> = q- x] (where q is a suitably-chosen 
wavevector) minimizes the energy. The third, forth and 
fifth terms vanish for the spiral state. Hence, the spiral 
state is a local minimum of the energy, as small perturba- 
tions around it would certainly raise the first two terms, 
and can only increase effect the remaining terms (as they 
are zero to begin with). The implication of this analy- 
sis is as follows: in the Double-Exchange Model there is 
a region of the zero-temperature (/i, Jh) phase diagram 
in which a spin spiral state is (at least locally) a sta- 
ble ground state. This state emerges on the high-/i (or, 
equivalently, 1ow-Jh) side of the continuous phase tran- 
sition line, on the other side of which the ferromagnetic 
state is the stable state. By using the fourth-order terms 
to restabilize the instability caused by the negative stiff- 
ness and including the effects of Jaf, we find that the 
wavevector a of the spiral is given by 



2 = 4J|(m 2 ~ J h + 32tt J h Jaf) 



1 - 16tt( J af / J h ) ' 
(2.14) 

where [j% = J H (1 — 327tJhJaf)- The approximate form 
holds for fi > fi c . 

Here, we note that if Jaf is greater than zero then the 
spin-spiral would not persist to arbitrary small carrier 
density, as the anti-ferromagnetic state becomes stable. 
Is such a continuous transition pre-empted by a first- 
order transition into a microphase-separated state? We 
explore this possibility in the next section. 
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FIG. 3: Arrangements of background spins on the 2x2 pla- 
quettes that are considered in the text. 



have meso-scale structure (magnetic and/or electronic) 
controlled by a competition between long-range inter- 
actions and interfacial energies. Prior work 0, has 
focused on the competition between the super-exchange 
and double-exchange coupling strengths, and has com- 
monly assumed the latter to be infinite, or at least very 
large. In the present setting, we are concerned with the 
entire range of double-exchange coupling strengths, but 
only with small super-exchange coupling strengths. 

In order to have the coexistence of the ferromagnetic 
and antiferromagnetic microphases associated with the 
microphase separation that we are considering, their 
thermodynamic and chemical potentials should coincide 
with one another These condition are necessary for 
microphase coexistence, but they are not sufficient in set- 
tings involving long-range interactions, such as those due 
to distinct charge-densities in the coexisting microphases, 
or interface energies associated with regions separating 
microphases. However, by examining the complement of 
the regions of the phase diagram that satisfy the afore- 
mentioned necessary conditions or are antiferromagnetic, 
we can locate the regions in which the homoegenous fer- 
romagnetic state or the textured state have a chance of 
being stable. (As we shall be limiting our consideration 
to the various types of antiferromagnetic ordering listed 
in Fig. |31 we may fail to exclude some regions of the 
phase diagram that we shall be calling ferromagnetic or 
textured.) 

We proceed by locating those regions of the (Jk,ii) 
phase diagram in which either the conditions for 
microphase-coexistance are satisfied or there is an an- 
tiferromagnetic state of lower energy than the ferromag- 
netic state. The single-electron band structure of the DE 
Hamiltonian (|2.1|l for ferromagnet (F), G-type (G) and 
A-type (A) antiferromagnet are given by For ferromag- 
net: 



III. PHASE SEPARATED AND CANTED 
STATES 

A. Collinear magnetic states 

In this section we compare the energies of sev- 
eral commonly-studied types of microphase-separated 
states jl6j in double-exchange magnetic systems, with the 
aim of comparing their stability relative to the spin spi- 
ral state. By microphase-separated we mean states that 



(F) ek = ±Jh ± 2tcosk x ± 2t cos k y — /i, (3.1a) 

(G) e k = ±\Jjl + 4t 2 (cos k x ± cos k y ) 2 - fi, (3.1b) 
(A) e k = ±2t cos k x ±J~J^+ 4i 2 cos 2 k y - fi. (3.1c) 

For both antiferromagnetic arrangements (G and A), all 
the eigenvalues are doubly degenerate. The energy den- 
sity and the electron density at zero temperature are 
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FIG. 4: Dependence of the ground state of the Double Ex- 
change Model on the carrier density and the Hund coupling 
at zero antiferromagnetic coupling. At low carrier densities 
and high Hund couplings ferromagnetism is favored. In cer- 
tain other regions, distinct homogeneous states are favored, 
including spin spiral (see Fig |SJ , A- and G-type antiferro- 
magnetism (see Fig. [3J. In yet other regions several types of 
FM/AFM microphase-separated states are favored. As ex- 
plained in the text, the boundaries indicated for microphase- 
separated states are in fact stability limits; the true bound- 
aries must lie within these stability limits. Only in undotted 
region can one be certain that the spin-spiral state is the fa- 
vored state. (Color version available online.) 

given by 

E = J 0e(e t ) e ,, (3.2a) 

BZ 

N = J 06(6,), (3.2b) 

BZ 

where for the antiferromagentic cases the Brillouin zone 
should be halved (in each direction), relative to the fcr- 
rmoagnetic case. 

From these ingredients we find numerically the lines 
in the (/i, ( J H )) plane at which the A-AFM/FM and 
G-AFM/FM phase transitions occur. On these phase 
boundaries, for each pair of competing states we calcu- 
late a pair of lines n c ( Jh) corresponding to the density 
of each state. On the (Jn,n) phase diagram these lines 
bound the regions depicted in Fig. ^inside which the two 
competing microphases have a chance of coexisting. 

B. Uniform canted magnetic states 

Next we consider the instability of the ferromagnetic 
state with respect to the canted ferromagnet state, in 
order to determine the phase boundary between them. 
Repeating the proceedure of the previous subsection, we 



FIG. 5: Background spin configuration in the canted ferro- 
magnetic state. 

find the single-electron band structure of Eq. I|2.1|l 

e/j = —2t cos k y — /i 

I (3-3) 

± V J h + 2<2 (! + cos2fc x ) ± 4J H i cos k x cos 6, 

where 9 is the canting angle of the background spin with 
respect to the 2— axis, as shown in Fig. [5] To find the 
phase transition line in the (//, Jh) diagram, we determine 
when the stiffness of the ferromagnetic state with respect 
to canting becomes zero. The resulting phase-boundary 
line is shown in Fig.0J From the diagram, it appears that 
the canted state is always pre-empted by a microphase- 
separated state. 

IV. CONCLUDING REMARKS 

To study smoothly varying textures of the Double Ex- 
change Model, we have followed the familiar program of 
expanding the free energy of the model in powers of gradi- 
ents of the background spin texture and tracing out the 
electronic degrees of freedom. This program has wide 
applicability to the study of long wavelength patterns in 
lattice models such as the Double Exchange Model. The 
main result of our paper is the phase diagram^ for the 
Double Exchange Model, which we have obtained analyt- 
ically as a function of the carrier density and the Hund 
coupling Jh between the carrier spins and the lattice of 
classical backgropund spins. Through the application of 
this program, we find that the spin-spiral state is in- 
deed a stable state for low carrier-densities and has a 
continuously- varying wavevector a ~ \fi — nd 1 ^ 2 - By di- 
rect diagonalization we also find that the transition from 
the ferromagnetic state to the canted state is essentially 
pre-empted by phase separation into different types of 
antiferromagnetic states. 
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APPENDIX A: SIMPLE MODEL OF SPIN 
SPIRAL TEXTURE IN TWO DIMENSIONS 



In this appendix, we solve for the exact ground-state 
energy of the DEM for the spin spiral state as the back- 
ground spin texture. We confirm that in two dimensions 
the critical line for the transition lies at \x — Jh, as stated 
in Sec. Ill Bl The integrating-out of the electronic de- 
grees of freedom, which we have carried out in the main 
text to determine the ground-state energy, becomes more 
transparent through this example, in which the process 
is nonperturbative. The simplification ensues when on 
restricts attention to a specific class of background spin 
configurations, viz., spin spirals, which have the form 




FIG. 6: Spin spiral configuration of background spins, with 
cone angle \ and wavevector a. 



0(x) = x ; (f>{x) 



where a is a constant wavevector. 
A(x) = a/2cosx, A.(x) — iasinx 



(Al) 



With this choice, 
and (<9n) 2 = 



a 2 sin x- The Hamiltonian 12.41 then reduces to 



H2D - 



(ipi(k) yt{k))( k 2 + a- fccosx + a 2 /4- J H - A* 

\ —a ■ k sinx k 2 



—a ■ k sin% 
fccosx + a 2 /4 + Jh 



fifth) 



4tt 2 



(A2) 



On diagonalizing this Hamiltonian and calculating the 
effective energy for the spin spiral state, we find that the 
critical chemical potential is equal to the Hund coupling. 
When ji is larger than Jh the ferromagnetic state is un- 
stable with respect to the formation of a spin spiral state. 
For x = tt/2 and fj, < Jh the energy has the form 



( Jh + M)2 I a 2 J h 



32ttJ h 



-a 



(4 -A + f 

256ttJ h 



For fi > Jh, minimizing E e s with respect to a gives 



4J h (m 2 -J h ) 



(A4) 



(A3) which determines the pitch of the stabilizing spin spiral. 



APPENDIX B: EVALUATION OF FEYNMAN DIAGRAMS 



In this appendix we present a sample diagram calculation by considering the two diagrams of lowest order in the 
gradient expansion (see Sec. Ill B|l . The Matsubara Green functions for aligned (tp) and anti-aligned (ip) electrons are 



G v (p,ip n ) 



1 



ipn - (p 2 - M - Jh) 



1 



ipn - (p 2 - M + Jn) 



(Bl) 
(B2) 



The two diagrams that contribute to the ferromagnetic stiffness are diagrams 1 and 2 in Fig. |SJ The amplitude 
corresponding to the first diagram is 



TjYl / dc l{G4,(q,iqn) + G v {q,iq n )}-{d/3n-df}h). 



(B3) 
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FIG. 7: Band structure in the presence (solid line) and ab- 
sence (dashed line) of the spin-spiral magnetic state in two 
dimensions. Notice that the Fermi surface just touches the 
bottom of the anti-aligned band for the optimal spin-spiral 
state. 



where dp stands for d p/ (2ir) d . The integral over the internal momentum results in the following expression, which 
is essentially the the sum of the volumes of two d— dimensional spheres of respective radii — Jh and + 7h : 

2 ~iw^~ { e(M + Jli)(/i + Jn)d/2 + e(/ ' - Jh)(m - Jn)rf/2 } • (B4) 

The amplitude corresponding to the second diagram is 



4^ / 3qakG v (q + k,t? n )G^(q,t« B )A /1 (k)A:(k)( g + */2) |t (« + */2) v (B5) 

= IV fa A„(k)A;(k)(g + fc/2) M ( g + fc/2)„ 

f - J q (iq n - (q 2 - (i - Jn))(iQn - (|k + q| 2 - M + Jk) ' 1 j 

On applying the standard Feynman trick and simplifying, the previous amplitude becomes 

I V fdz ^ (k)A - (k)(g + fc/2) - (g + fc/2) - f (B7) 

[(l-z)( ig „-( g 2 -M- l /H)) + ^^n-(|k + q| 2 - Ai + J H )] 2 



A„(k)A*(k)^ g!/ 



^ij^ "'o [i<7„-<7 2 +M-£: 2 z + /c 2 z 2 + Jh(1-2z 



]2 ' 



(B8) 



where in the final step we have dropped terms of higher order in k, as they do not contribute to the stiffness. 
This follows from the obsevation that A* A Q = dpii ■ d^h. Reversing the order of summation and integration, and 
simplifying further by noting that the denominator sums to a <5-function in the zero-temperature limit, we obtain 

J dqdk J dzA tl (k)A* L ,(k)(q tl q„)S[-q 2 + [i-k 2 z + k 2 z 2 + Ml-2z)] (B9) 
~9[ M - JgK/x - J H ) 1+d/2 + 9[J H + m](Jh + /i) 1+rf / 2 f 

= w/ 2 j H r[2 + d/2] J dx • « ' (B10) 

which simplifies to the expressions l|2.10|l in the text. The other diagrams, which are more complicated are nevertheless 
evaluated in a similar fashion, and since the technique is rather well established and used commonly, we find it 
redundant to present the details here. 
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